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ABSTRACT 

Context. INTEGRAL, the European Space Agency's 7-ray observatory, tripled the number of super-giant high-mass 
X-ray binaries (sgHMXB) known in the Galaxy by revealing absorbed and fast transient (SFXT) systems. 
Aims. In these sources, quantitative constraints on the wind clumping of the massive stars could be obtained from the 
study of the hard X-ray variability of the compact accreting object. 

Methods. Hard X-ray flares and quiescent emission of SFXT systems have been characterized and used to derive wind 
clump parameters. 

Results. A large fraction of the hard X-ray emission is emitted in the form of flares with a typical duration of 3 ks, 
frequency of 7 days and luminosity of 10 36 erg/s. Such flares are most probably emitted by the interaction of a compact 
object orbiting at ~ 10 R* with wind clumps (10 22-23 g) representing a large fraction of the stellar mass-loss rate. The 
density ratio between the clumps and the inter-clump medium is 10 2-4 in SFXT systems. 

Conclusions. The parameters of the clumps and of the inter-clump medium, derived from the SFXT flaring behavior, 
are in good agreement with macro-clumping scenario and line driven instability simulations. SFXT have probably a 
larger orbital radius than classical sgHMXB. 
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1. Introduction 

Stellar winds have profound implications for the evolution 
of massive stars, on the chemical evolution of the Universe 
and as source of energy and momentum in the interstellar 
medium. 

Photons emitted by massive stars directly transfer mo- 
mentum to the stellar wind through absorptio n in numerous 
Dopp ler shifted optically thick spectral lines (jCastor et alj 
Il975l ) and drive the wind to highly supersoni c speeds. Such 



line driven stellar win ds are very unstable (jOwocki et alJ 
1988; [Feldmeierl Il995) . The collisions between high speed 
and low density material with slower gas trigger strong 
shocks, form high density contrasts and wind clumps and 
lead to thermal X-ray emission. 

There are multiple observational lines of evidence for 
wind clumping in massive stars from w ind line profiles 
(|Bouret et a lJ l2005t iFullerton et al.ll2006lh d iscrete absorp- 
tion co mponents (Prima fe Howarthl 19861). variable line 
profiles dLepine fe MoffatJll999t iMarkova et all 12005ft. po- 
larimetry (ILupie fc Nordsieckl[l9 87: Dav ies et al.ll2007l ). X- 
rays continuum dCassinelli fc Qlsonlll979l ) and line emission 
(jOskinova et al.ll2006f) . ~ 

Wind clumping has an important effect on mass-loss 
rate diagnostics that depends on the square of the wind 
density. It leads to a reduction of the stellar mass-loss rate 
estimates by a factor fy ' 5 wher e fy is the clump volu me 
filling factor (| Abbott et al.lll98ll : iFulTerton et al.N2"006h . A 
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reduction by a fact or > 3 becom es problematic for mas- 
sive star evolution (jHirschil 120071 ) . Optically thick clumps 
may however help to reconcile wind clumping and usual 
mass-loss rates by reducing sp ectral features and free-free 
emission (|Oskinova et al.ll2007r ). 

Indirect measures of the structure of massive star winds 
are possible in binary systems through the analysis of the 
interaction between the stellar wind and th e companion 
or its stell ar wind. In col l iding wind binaries, ISchild et al.l 
(|2004l ) and iPollock et ail (|2005f ) have pointed out that the 
observed X-ray column densities are much lower than ex- 
pected from smooth winds an d that clump ed winds are 
the probable solution (see also iPittardl 12007). In wind-fed 
accreting High-Mass X-Ray Binaries (HMXB), clumping 
has been invoked to explain the orbital variabilit y of X-ray 
line profiles emitted by the pho to-ionized wind (jSako et al.l 
1200a Ivan der Meer et alll2005h . 

In this paper we aim at studying the clumping of stel- 
lar wind in HMXB using the hard X-ray variability ob- 
served by the IBIS / ISGRI instrument (jUbertini et aljfeooa 
iLebrun et al.ll2003l) on board the INTErnat ional Gamma- 
Ray Astrophysics Laboratory (INTEGRAL, IWinkler et al.l 
l2003h which could be used to probe clump parameters and 
density con trasts in the stellar wind. This study follows a 
first paper (|Levder et al.ll2007h interpreting the flaring be- 
havior of IGR J08408— 4503 in term of wind clumping. 

Classical wind-fed super-giant HMXB (sgHMXB) are 
made of a compact object orbiting within few stellar radii 
from a super-giant companion. Recently INTEGRAL al- 



most tripled the number of sgHMXB systems known in the 
Galaxy and revealed a much more complex picture with two 
additional fa milies of sources: the hig hly absorbed persis- 
tent systems ([Walter et al.112003 . |2006| ) a nd the super-giant 
fast X -ray transient (SFXT) systems (jNegueruela et al.l 
l2006bh . 

The highly absorbed systems have orbital and spin 
periods similar to those of classical Roche-lobe under- 
flow sgHMXB, however the absorbing column densities are 
much higher than o bserved on average in classical systems 
(| Walter et al.l [200^ 1. The fast transient systems are char- 
acterized by fast outbursts, by a low quiescent luminos- 
ity and by super-giant OB companions (jSguera et al.i r2006: 
iNegueruela et al.ll2007l ). 

The sample of sources and the INTEGRAL data analy- 
sis are described in Sect. 2. Their hard X-ray variability is 
discussed in the context of clumpy stellar winds in Sect. 3. 
Finally Sect. 4 summarizes our principal conclusions. 



2. Sources and data analysis 

2.1. Source sample 

Many sources have now been proposed as candidate super- 
giant fast X-ray transient based on their hard X-ray vari- 
ability characteristics, and, for a subset of them, opti- 
cal counterpart spectral type. Contrasting statements have 
however been made on specific sources for what concerns 
their persistent or transient nature. In the frame of the 
current paper we have considered all SFXT candidates to- 
gether with several persistent and absorbed super-giant 
HMXB for comparison. 

Among them, we specifically excluded known Be 
systems (IGR J11305-6246, AX J1700. 2-5129), sources 
detected only once by INTEGRAL (IGR J16283- 
4843, IGR J16358-4726), blended INTEGRAL sources 
(IGR J17407 -2808), the long period systems IGR 
J11215-5952 (ISidoli et al.l l2006li and AX J1749. 1-2733 



(jZurita Heras k Chatvll2007h an d the sgB[e] system IGR 
J16318-4848 (fWalter et al.ll2006l) . 

We analyzed the available INTEGRAL data for 12 can- 
didate SFXT that have large variability factors (table [T]) 
and compared them with classical and absorbed sgHMXB 
system that have typical variability factor < 20. 

The 12 classical sgHMXB systems (4U 1700-37, 4U 
1538-52, SAX J1802. 7-2017, Cyg X-l, XTE J1855-026, Vela 
X-l, 2S 014+65, IE 1145.1-6141, GX 301-2, 4U0900-40, 
4U1907+09, XTE J1855-026) have orbital radii ranging be- 
tween 1.5 and 2.5 stellar radii. 

The 7 abs orbed sgHMX B system s (IGR 
J16320-4751 (iRodriguez et all I2006D . IGR 
J16393-4641 (iBodaghee et al.l 12006b. IGR J1 7252-3616 
(IZurita Her as et alj l2006f), IGR J18027-2016 dHill et all 
I2005h. IGR J1 9140+0951 ([Rodriguez et all |2005j), IGR 
J16418-4532 ([Walter et all l2006h. IGR ,100370+6122 
(jin't Zand et al.ll2007l )f have orbital periods < 15 days and 
probable super-giant companions i.e. small orbital radii as 
well. 



2.2. INTEGRAL data and analysis 

The sources of the sample are located along the galactic 
plane that has been heavily observed by INTEGRAL. All 
public data available in March 2007 are considered in this 



study. The INTEGRAL imager IBIS has a large field of 
view (FOV) of 29° square with a good spatial resolution 
of 12'. Several thousands pointings including at least one 
source of the sample in its FOV and with an exposure longer 
than 600 s were selected, spanning times between Jan. 11, 
2003 (revolution 30, MJD 52650) and Dec. 2, 2005 (revolu- 
tion 383, MJD 53709). The total elapsed exposure time on 
each source goes from 52 to 127 days. The observations are 
not equally distributed in time due to scheduling reason. 

The ISGRI data are reduced using the INTEGRAL 
Offline Scientific Analysis (OSA) version 6.0 software pub- 
licly released by the INTEGRAL Science Data Centre 
(jCourvoisier et all 120031 ) . Individual sky images for each 
pointing have been produced in the energy band 22-50 keV. 
The detection of the sources of the sample is forced in each 
individual images and the source count rate extracted. 

Sky maps gathering all pointings within the same rev- 
olution are also built. To build mosaic images with longer 
exposure time, each source is treated separately. Sky mo- 
saics are built combining consecutive pointings in which 
the source is not detected at least at a 5.1c level. In several 
cases, when the source is not detected in several revolu- 
tions that are close in time, they are also mixed together to 
create one deep-exposure mosaic. Finally, one last mosaic 
is created combining all pointings where the source is not 
detected in order to constrain the lowest hard X-ray emis- 
sion level. The light curves are built from the sky images 
considering only the pointings where the source is located 
less than 14° away from the FOV centre. 

The source count rate is extracted with the tool 
mosaic_spec (version 1.4) that is part of the OSA pack- 
age. Detection of the source in mosaics are considered at 
a 6cr level or higher. Therefore, each source light curve is 
built mixing short- and long-term activity. The 22-50 keV 
count rates can be converted into physical fluxes with the 
relation 1 Crab = 117ct/s ~9x 10~ 9 ergcm~ 2 s _1 . 

Because of the INTEGRAL dithering operations, the 
source could move out of the FOV and returns a few point- 
ings later and still be detected. Source flares have therefore 
been detected by requiring a minimum of 25 ks of inactivity 
be tween them. Examp les o f flare lightcurves a r e pres ented 
in lSguera et~ai1 (|2006l ) and lZurita Heras et all (|2007t) . 

Flare duration of the order of a single INTEGRAL 
pointing (2 ks) have been observed in all sources, (excepting 
IGR J16465— 4507, but see below). Their typical duration 
is 3 ks. 

Fewer longer (> 15 ks) flares have also been detected. As 
many flares are detected close to the sensitivity limit of the 
instrument, the flare duration is often difficult to evaluate 
and long flare could often also be considered as a series of 
shorter flares. We inspected each long flare carefully and a 
number of cases were found: 

- In several sources (IGR J17544-2619, XTE J1739- 
302, SAX J1818.6-1703) some long flares were found 
to be composed of very strong short flares surrounded 
by activity at much lower level. Such flares were finally 
accounted in table Q] as one short flare. Interestingly 
the detailed analysi s of a flare of IGR J08408-4503 
(jLevder et alJ 120071 ) also revealed that it was sur- 
rounded by emission at a lower level during several 
hours. This could be a common behavior, only detected 
by INTEGRAL for the brightest flares. 



Table 1. List of SFXT candidates with quiescent flux F q , flaring characteristics: range of flare count rates Ffi, number of 
short (< 15 ks) and long (> 15 ks) flares Nfi, the average and range of durations for short flares and the list of durations 
for long flares tfi and source observing elapsed time T b s . 



Source 


Spectral 


Distance 


F q 
ct/s 


Ffi 


Nfi 


tfi [short] 


tfi [long] 






Type 


kpc 


ct/s 


short+long 


ks 


ks 


days 



SFXT systems 



IGRJ08408-4503 


08.5 lb 1 


2.8 1 


< 


1 


2.1-3.9 


2+0 


3.6 








52.0 


IGRJ17544-2619 


09 lb 2 


2-4 2 





06 


4.2-24 


8+0 


2.5 


(2 


-4.3) 




127.0 


XTEJ1739-302 


08.5 Iab(f) 3 
09-Bl I 4 


1.8-2. 9 3 





08 


3.0-28 


12+1 


4.2 


(2 


-8) 


50 


126.4 


SAXJ1818.6-1703 







18 


5.2-45 


11+0 


2.9 


(2 


-6) 




76.9 


IGRJ16479-4514 


OB I s 







2 


2.5-19 


27+11 


3.6 


(2 


-14) 


(16-35) + 84 


67.0 


AX J1841. 0-0536 


BO I 6 




< 


1 


3.7-15 


4+0 


5.8 


(2 


-13.1) 




51.9 


AXJ1820.5-1434 


B 4 




< 


1 


3.4-5.3 


4+0 


3.9 


(2 


-9.6) 




59.4 



Intermediate systems 



IGRJ16465- 


-4507 


B0.5 I 7 




0.1 


2.5-6.9 


0+3 








19, 25, 45 


66.7 


AXJ1845.0- 


-0433 


09 la 8 


3.6* 


0.2 


4.0-6.2 


6+0 


4.0 


(2 


-14.3) 




55.2 


IGRJ16195- 


-4945 


OB I 9 


7 9 


0.2 


2.8-4.8 


6+0 


2.2 


(2 


-3.3) 




71.8 


IGRJ16207- 


-5129 


B0 I 4 


4 4 


0.4 


2.8-9.2 


9+2 


4.3 


(2 


-11) 


18, 25 


73.7 


XTEJ1743- 


-363 






0.5 


4.1-9.2 


16+3 


2.5 


(2 


-6.6) 


21, 45, 61 


122.9 



Ref.: Uevder et al.lil2007h. iPellizza et all (120061 ). jlNeeueruela et all (|2006ah. INegueruela fc Schurchl ([20071 ). IChatv fc al.llj2007n . 
^ iNespoli et al.l i|2007l ), INegueruela et al l (|2005l ). ICoe et al.1 (|1996l ). ISidoli et al.1 (I2005T T 



- In XTE J1739-302 and SAX J1818.6-1703 a few long 
flares were clearly the result of two short flares separated 
by less than 25 ks. Such flares were accounted as two 
short flares in table [TJ 

- In IGR J16207-5129, XTE 1743-363 and IGR 
J16465— 4507, long flares were found to be sequences 
of detection in individual pointings, at the limit of de- 
tectability, separated by pointings where the source was 
in the FOV but too weak to be detected. When de- 
tected, the source count rate was similar to the sen- 
sitivity limit at the source off-axis angle. We conclude 
that those flares are particularly long periods of activity 
at the limit of the instrumental sensitivity. Such cases 
were kept in the "long flare" category in table [TJ These 
sources are likely classical sgHMXB, with average fluxes 
just below the IBIS sensitivity. Such long activity peri- 
ods cannot be considered as real flares. 

- IGR J16479-4514 is very active with one flare detected 
per day on average. 11 long flares were detected. 8 of 
them feature few very significant detections at pointing 
(2 ks) level separated by periods of 10 — 20 ks during 
which the source was not detected. These flares could in 
fact be split in shorter flares with duration < 8 ks. Three 
long flares are detected close to the detection limit and 
could also be represented as sequences of shorter flares. 
The large number of long flares detected in this source 
is therefore probably spurious and related to the high 
flare frequency. As a detailed analysis of the numerous 
flares of this source is beyond the scope of this paper 
we kept these long flares separated in table [TJ 



2.3. Variability and Source Classification 

Table [TJ lists the sources together with their quiescent count 
rate (F q ), typical flare count rate {Ffi), number of short 
(< 15 ks) and long (> 15 ks) flares (Nfi), range of flare 
durations (tfi) for short and long flares and total source 
observing elapsed time (T b s ). The sources are ordered from 
high to low variability factor (Ffi/F q ). 

The source observing time T t, s is the sum of the elapsed 
time of all pontings with the source within 14° of the FOV 
center. As the instrument effective area decreases between 
the borders of the fully and partially coded fields of view, 
the probability to detect a flare effectively decreases when 
the source gets outside of the fully coded field of view. The 
effective observing time for flare detection can be estimated 
as 0.6 T i, s . The effective time period between two flares is 
typically T — 7 days on average but varies between 1 day 
in IGR J16479-4514 and two weeks in IGR J08408-4503 
(this is uncertain as only two flares were detected). 

The flare count rates (Ffi) are averaged over the dura- 
tion of each flare. Peak count rates could be significantly 
larger than listed in table [TJ espec ially for flares sho rter 
than the pointing duration (see e.g. lLevder et al.|[2007f ). 

The sources have been separated in two categories. The 
SFXT includes systems featuring hard X-ray variability 
by a factor > 100. "Intermediate" systems are candidate 
SFXT with smaller variability factors that could be com- 
pared with those of classical systems. The dividing line be- 
tween the two source categories is not very well defined. 
From the variability point of view, sources closer to the 
bottom of the table are more similar to classical sgHMXB. 



We only have one really long flare remaining. A 50 ks 
flare has been detected in XTE J1739— 302 as a continuous 
period of activity. It is very structured and made of at least 
3 different peaks. It could be interpreted as a close sequence 
of shorter flares. 



3. Discussion 

3.1. Hard X-ray luminosities 

The distance to the SFXT systems has been evaluated in 
a few cases (jLevder et al.l 120071 : INegueruela et al.l l2006at 



iPellizza et alll2006l ). They range between 2 to 7 kpc with 
large uncertainties. We will assume, for the rest of the dis- 
cussion a distance of 3 kpc. The average count rate ob- 
served during flares lies between 3 and 60 ct/s which trans- 
lates to hard X-ray luminosities of (0.2 — 4) x 10 36 erg/s. 
Such luminosities are not exceptional for sgHMXB but 
very significantly larger than the typical X-ray luminos- 
ity of single massive s tars of 10 30-33 erg/s at soft X-rays 
(jCassinelli et al.lll98ll ). The observed non thermal hard X- 
ray emission is therefore characteristic of accretion on a 
putative compact object and is, for now, the main reason 
to assume that SFXT are binary systems. 

As the sources are flaring at most once per day, their 
average hard X-ray luminosity is very low, reaching (0.2 — 
4) x 10 34 erg/s. It is therefore very unlikely that those sys- 
tems have average orbital radius lower than 10 13 cm i.e. 
~ 10 i?*. One expects orbital periods larger than 15 days 
and underflow Roche lobe systems (note that no orbital 
period has yet been derived in any of these systems). 

The average hard X-ray luminosity of the SFXT systems 
in quiescence is < 6.6 x 10 33 erg/s (which corresponds to 0.2 
ct/s). This is an upper limit as the mosaics used to measure 
those quiescent fluxes most probably contain faint flares, 
not detected during single pointings. The INTEGRAL data 
alone do not exclude that there is no quiescent hard X-ray 
emission in these systems. 

3.2. Wind clumps 

The interaction of a compact object with a dense 
clump formed in the wind of a massive companion leads 
to increased accre tion rate and hard X-ray emission 
(jLevder et al.ll2007t) . 

The free-fall time from the accretion radius R a = 
2 x 10 10 cm towards the compact object is of the order 
of (2 — 3) x 10 2 s. As the intrinsic angular momentum o f 
the accreted gas is small (jlllarionov fc Beloborodovll200l[ ) 
the infall is mostly radial (down to the Compton radius) 
and proceeds at the Bondi-Hoyle accretion rate. 

The accretion could slow down if the wind clumps have 
internal turbul ence or harbor significant intrinsic angu - 
lar momentum (|Theuns et al.| [l996: Krumho lz et al.l l2005). 
This is however unlikely in a highly supersonic wind, and 
supported by the very sharp X-ray flare cutoff observed on 
time scale of the order of few 100 s in some of these sys- 
tems (jZurita Heras fc Waited 120071 ; iGonzalez-Riestra et al.l 
1200411 . 

With a duration of tfi = 2 — 10 ks, the observed short 
hard X-ray flares are significantly longer than the free-fall 
time. The flare duration is therefore very probably linked 
with the thickness of the clumps which, for a clump radial 
velocity V c i = 10 8 cm/s, is h d = V c i x tfi ~ (2 — 10) x 
10 11 cm. 

The average hard X-ray luminosity resulting from an 
interaction between the compact object and the clump can 
be evaluated as Lx — e M acc c 2 /tfi (where e ~ 0.1) and the 
mass of a clump can then be estimated as 

M cl = {Rci/Raf M acc = (R c i/R a ) 2 L x t f i/(e c 2 ) 

where R c i is the radius of the clump perpendicular to the 
radial distance. In the case of a spherical clump, 



M r , = 



( L 



If N is the rate of clumps emitted by the star, 
the observed hard X-ray flare rate is given by T^ 1 = 
N (R 2 cl / AR 2 rb ). The rate of mass-loss in the form of wind 
clumps can then be estimated as 



M rl = MrlN = 



/10d 



4R 2 orb L x tfi 
R 2 T e c 2 

L x t f i 



V T 10 36 erg/s3ks, 

For a (3 = 1 velocity law and spherical clumps, the num- 
ber of clumps located between 1.05i?* and R or b can be eval- 
uated as 



N = N (t(R orb )-t(1.05R*)) 

10 d\ /3 ks\ / R or b 



3 ksY 



3.8 x 10 3 , 



7 x 10~ 14 g cm" 3 



x 



\ 10 36 erg/s 



llL 
3 ks 



7.5 x 10 21 g. 



T J \ t f i J \10 13 cm 

where t(r) is the wind flight time (jHamann et al.ll200lh . 
Assuming spherical clumps, the clump density at the 

orbital radius is p d = (ipirf^T^ 
and the corresponding homogeneous wind density is 
Ph = M c/ /(4tt R 2 orb V a ) = (^To^iTi^) 1-5 * 
10~ 15 g cm~ 3 . The clump volume filling factor at the 
orbital radius is f v = = ( W^ ) 0.02 and 

the corresponding porosity length (jOwocki &: Cohen |[2006; 
lOskinova et aLll2007ft is h = = (y^) 15 x 10 12 cm. 

If the density of a clump decreases with radius as 
r~ 2/3 and its mass remains constant, the averaged homoge- 
neous wind density within R b s is 7j/T = NM c i/(^nR 3 rb ) = 

10 T d io^?rg/ s rfe) 7 x 10_15 S cm ~ 3 and the average 
clump volume filling factor and porosity length could be 
estimated as 0.1 and 3 x 10 12 cm, respectively. 

The variety of tfi, T and Ffi that are observed proba- 
bly reflects a range of clump parameters and orbital radii. 
Several of the average clump parameters estimated above, 
in particular the clump density, filling factor and porosity 
length do not depend on the orbital radius, which is un- 
known, and only slowly depend on the observed quantities. 

These average par ameters match t he ma cro-clumping 
scenario proposed by lOskinova et al.l (|2007f) to reconcile 
clumping and mass-loss rates. Their model depends on a 
free parameter Lq = L(r)(r 2 V(r)/V ^(oo))" 1 / 3 , where L(r) 
is the clump separation in unit of Our average clump- 
ing parameters correspond to L = 0.35. The number of 
clumps derived above is also comparable to evaluat i ons b y 
iLenine & MoffaH (|1999f ) and lOskinova et all (I2006L 120071 ). 
The volume filling factor and the c lump mass-lo s s rate 
are also similar to those derived by iBouret et al.l |2005) 
from the study of ultraviolet and optical line profiles in two 
super-giant stars. 

The column density through a clump can also be esti- 
mated as N H = -S^- = ( , n £ x i Ur 

clumps remain optically thin in the X-rays. 

The activity surrounding some of the most significant 
flares may be related to tidal effects on the clumps them- 
self and induced turbulence. The long flare that has been 
observed in XTE J1739— 302 is structured with many sub- 
peaks and may correspond to the clustering of smaller wind 
clumps perhaps related to the dissolution of pancake or 
shell shaped wind clumps. 



5 x 10 22 cm- 2 . The 



3.3. Missed flares 

The mass spectrum of the wind clumps is completely un- 
known. Only the brightest hard X-ray flares in SFXT 
and "intermediate" objects could have been detected by 
INTEGRAL. Bright and short as well as long and faint 
flares could have remained undetected. 

By definition, the contribution of missed flare to the 
average source count rate could not be differentiated from 
the quiescent emission. The average contribution of unde- 
tected flare, which is not larger than the quiescent count 
rate listed in table [TJ can be compared to the contribution 
of the detected flares. 

In two SFXT (IGR J16479-4514, AX J1820.5-I434) 
the average count rate is more than 3 times larger than 
the quiescent count rate. In these cases, the global contri- 
bution of undetected flares is at most a third of the con- 
tribution of the detected flares. In three additional sources 
(IGR J16465-4507, XTE J1739-302, AXJ 1841.0-0536) 
the average count rate is slightly larger than the quiescent 
count rate and the global contribution of undetected flares 
is at most similar to the contribution of the detected flares. 
In these 5 sources, INTEGRAL thus detected a significant 
fraction of the accretion luminosity. 

In the two remaining SFXT (IGR J17544-2619, SAX 
J1818.6— 1703) and in the other "intermediate" sources, the 
average count rate is compatible with the quiescent count 
rate. One cannot therefore exclude that the contribution 
of undetected flares to the average hard X-ray luminosity 
could be larger than that of detected flares. 

A large fraction of the hard X-ray emission detected in 
SFXT is therefore generated by the accretion of the clumps 
described in Sect. 13.21 

A population of small clumps (with tfi < 300 s) could 
have remained almost undetected if their interaction lumi- 
nosity is not larger than 10 36 erg/s. Their mass would be 
of the order of 10 19 g. Short flares would need to be very 
frequent (T < 1 day) for those clumps to contribute signifi- 
cantly to the mass-loss rate and should be easily detectable 
in X-ray observations of a sample of sources. 

Undetected long (tfi > 50 ks) and faint (Lx < 
10 35 erg/s) flares would be the signature of extremely large 
(R c i > i?*) and massive (10 24 g) clumps with a density of 
the order of the equivalent homogeneous wind density. Such 
structures would be interpreted as wind in-homogeneities 
with low mass-loss rate rather than by clumping. 

3.4. Inter-clump medium 

The variation of the observed X-ray flux between flares and 
quiescence provides in principle a direct measure of the den- 
sity constrast between the wind clumps and the inter-clump 
medium. As part of the quiescent count rate could be re- 
lated to undetected flares or to other emission, these density 
contrasts are lower limits. 

As the hard X-ray flux variability observed in SFXT 
(table [T|) reaches factors > 100 when a quiescent count 
rate could be firmly detected with INTEGRAL, the density 
contrast is at least of a factor 100. In the "intermediate" 
sources INTEGRAL observed variability factors larger than 
20. Additional constraints are available for the sources ob- 
served at higher sensitivity in the X-rays: 

- In IGR J08408-4503, iLevder et all (|2007t ) have anal- 
ysed SWIFT XRT observations which provided a very 



faint quiescent luminosity of 2 x 10 32 erg/s in the form 
of a soft, probably thermal, spectrum. In this source the 
wind density contrast was larger than 10 4 . 

— The flaring activit y of IGR J17544— 2619 was o bserved 
by XMM/Newton ([Gonzalez- Riestra et al]|2004h . These 
flares were short and very structured on time scale of 
100 s. The variation of the count rates observed between 
quiescence and average or peak flare are of a factor 20 to 
40. Similar and stronger variations have been observed 
by INTEGRAL ([Walter et al.ll2006h . 

- IGR J16479-4514 was observed by XMM/Newton and 
it was found 20 times f ainter than the ave rage flux ob- 
served by INTEGRAL ([Walter et al.ll2006h . The overall 
variability observed in this source is of the order of 100. 

— In AX J1845. 0—0433, the flare and quiescent fluxes 
were observed simultaneously by XMM-Newton. The 
quiescent spectrum is clearly dominated by accretion 
(|Zurita Heras fc Walteill2007li . A density contrast of 15- 
25 was detected. 

- IGR J16465- 4 507 w as observed by XMM-Newton 
([Walter et alJ l2006h and IGR J1619 5-4945, 
IGR ,116207 -5129 (iTomsick et alJ 12006ft . AX 
J1841.0-0536 (lHalpern et al.ll2004D . SAX J1818-1703 
(lin't Zand et al.ll2006D and XTE 1739-302 ([Smith et alJ 
2006) were observed by Chandra. In all cases, the X-ray 
flux is similar to the average quiescent fluxes observed 
by INTEGRAL (table Q]) . 

Density contrasts of > 10 2 ~ 4 and 15-50 have been 
observed in SFXT and "Intermediate" sources, respec- 
tively. The density contrast is larger in SFXT than in 
"Intermediate" and, of course, classical systems. Density 
contrasts are probably stronger when clumping is very ef- 
fective. 

Numerical simulations of the line driven instability 
([Runacres fc Owockl [2005) predict density contrasts as 
large as 10 3-5 in the wind up to large radii. At a distance 
of 10 i?*, the simulated density can vary between 10~ 18 
and 10 -13 g cm -3 and the separation between the density 
peaks are of the order of i?* . These characteristics are com- 
parable to the values we have derived from the high-energy 
observations. 

3.5. Where do classical sgHMXB fit in this picture ? 

Classical sgHMXB are characterized by small orbital radii 
Rorb = (1.5 — 2.5) i?*, and by flux variability of a factor 
< 20. Such variabilities were modelled in term of wind in- 
homogeneities largely triggered by the hydrodynamic and 
photo-ionisation effects of the accreting ob ject on the com- 
panio n and inner stellar wind (|Blondin et al.lll99ltlBlondin! 
Il994[ ). At small orbital radii, the companion is close to fill 
its Roche lobe, which triggers tidal streams. In addition the 
X-ray source ionizes the wind acceleration zone, prevents 
wind acceleration and generates slower velocities, denser 
winds, larger accretion radius and finally larger X-ray lu- 
minosities. Whether or not the stellar wind is intrinsically 
clumpy at low radius, the effect of the compact object on 
the wind is expected to be important. 

The main difference between SFXT and classical sgH- 
MXB could therefore be their orbital radius. At very low 
orbital radius (< 1.5 i?») tidal accretion will take place 
through an accretion disk and the system will soon evolve 
to a common envelope stage. At low orbital radius (~ 2 i?„) 



the wind will be perturbed in any case and efficient wind 
accretion will lead to copious and persistent X-ray emis- 
sion (10 36-37 erg/s). At larger orbital radius (~ 10 R*) 
and if the wind is clumpy, the SFXT behavior is expected 
as described above. If the wind clumps do not form for any 
reason, the average accretion rate will remain too low and 
the sources will remain mostly undetected by the current 
hard X-ray surveying instrumentation. 

4. Conclusions 

INTEGRAL tripled the number of super-giant HMXB sys- 
tems known in the Galaxy and revealed two new popula- 
tions: the absorbed and the fast transient (SFXT) systems. 
The typical hard X-ray variability factor is < 20 in clas- 
sical and absorbed systems and > 100 in SFXT. We have 
also identified some "intermediate" systems with smaller 
variability factors that could be either SFXT or classical 
systems. 

The SFXT behavior is best explained by the interaction 
betwe en the accreting compact object and a clumpy stellar 
wind (jin't Zand! 120051 : iLevder et al.ll2007th Using the hard 
X-ray variability observed by INTEGRAL in a sample of 
SFXT we have derived typical wind clump parameters. The 
compact object orbital radius are probably relatively large 
(10 i?*) and the clumps which generate most of the hard 
X-ray emission have a size of a few tenth of i?„ . The clump 
mass is of the order or io 22-23 g (for a column density 
of 10 22 ~ 23 cm -2 ) and the corresponding mass-loss rate is 
10 _< - 5-6 ^ M Q /y. At the orbital radius, the clump separation 
is of the order of i?* and their volume filling factor is 0.02. 
Depending how the clump density varies with radius, the 
average volume filling factor could be as large as 0.1. 

These parameters are in goo d agreement with the m acro 
clumping scenario proposed bv lOskinova et alJ ((2007). 

The observed ratio between the flare and quiescent 
count rates indicate density ratios between the clumps and 
the inter-clump medium which vary between 15 to 50 in 
"Intermediate" systems and 10 2-4 in SFXT. Such ratios 
and the observed clump densities are in reasonable agree- 
ment with the predictions of line driven instabilities at large 
radii (jRunacres & Owockil[2005h . 

The main difference between classical sgHMXB and 
SFXT could be the orbital radius of the compact object. At 
small orbital radius (R or b ~ 2 _R„) the systems are persis- 
tent and luminous. At larger radius and if wind clumping 
takes place the fast transient SFXT behavior is observed. 
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